Phase transition of clock models on hyperbolic lattice studied by corner transfer 

matrix renormalization group method 
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Two-dimensional lerromagnetic Ai'-state clock models are studied on a hyperbolic lattice repre- 
sented by tessellation of pentagons. The lattice lies on the hyperbolic plane with a constant negative 
scalar curvature. We observe the spontaneous magnetization, the internal energy, and the specific 
heat at the center of sufficiently large systems, where the fixed boundary conditions are imposed, 
for the cases A'^ > 3 up to A'^ = 30. The model with = 3, which is equivalent to the 3-state Potts 
model on the hyperbolic lattice, exhibits the first order phase transition. A mean-field like phase 
transition of the second order is observed for the cases A'^ > 4. When A'' > 5 we observe the Schottky 
type specific heat below the transition temperature, where its peak hight at low temperatures scales 
as A'^"'^. From these facts we conclude that the phase transition of classical XY- model deep inside 
the hyperbolic lattices is not of the Berezinskii-Kosterlitz-Thouless type. 

PACS numbers: 05.50.-|-q, 05.70.Jk, 64.60.F-, 75.10.Hk 



I. INTRODUCTION 

Two-dimensional (2D) lattice models with continuous 
local spin symmetry, such as the classical XY-model and 
the classical Heisenberg model on the square lattice, do 
not have finite magnetization when temperature is finite. 
This fact proved by Mermin and Wagner (Ij does not ex- 
clude the presence of phase transition of the Berezinskii- 
Kosterlitz-Thouless (BKT) type These well-known 

facts are based on analysis in the flat 2D plane. 

Quite recently, Baek et al. studied the XY model on 
the heptagonal lattice Q , which is one of the hyperbolic 
lattices constructed as a tessellation of heptagons on the 
hyperbolic plane, i.e., the 2D space with a constant nega- 
tive curvature @ . By way of the Monte Carlo (MC) sim- 
ulations for open boundary systems, they concluded the 
absence of phase transition, including that of the BKT 
type. Their result is in accordance with the thermody- 
namic property of the Ising model on the hyperbolic lat- 
tice, where there is no singularity in the specific heat as 
shown by d'Auriac et al. [6|. These observations on the 
hyperbolic lattice can be explained by the non-negligible 
effect of the system boundary 0, Q , which always has a 
finite portion of the system regardless of the system size. 

It should be noted, as pointed by d'Auriac et al, that 
the presence of the ordered phase is not excluded in the 
region far from the boundary Q, although the area of 
such an ordered region is negligibly small compared with 
the whole system on the hyperbolic lattice. The situ- 
ation is similar to that of the statistical models on the 
Cayley tree, where its deep inside can be regarded as the 
Bethe lattice @. Shima et al. studied the Ising model 
on the hyperbolic lattice by the MC simulations, and ob- 
served the mean-field like phase transition deep inside 
the system [13, [HI- The mean field behavior is in ac- 
cordance with theoretical studies of phase transition in 
the infinitely large hyperbolic lattices [l^, [l^ . It can be 



expected that such an order also appears in the case of 
the XY-model and the clock models. 

In this paper we study A^(> 3)-state clock models on 
the pentagonal lattic e [Til u p to iV = 30 by use of the 
CTMRG method [lljlajp] modified for systems on the 
hyperbolic lattices IsTlot. The internal energy and the 
spontaneous magnetization at the center of sufficiently 
large systems are calculated numerically. In order to 
judge the presence of an ordered state deep inside the 
system, we impose the ferromagnetic boundary condi- 
tions at the beginning of the iterative calculation of the 
CTMRG method. As we show in the following, the ob- 
tained results support the existence of the mean-field like 
phase transition for all the A'^ even in the limit A*" — > oo, 
where the system coincides with the classical XY model. 

In the next section we introduce geometry of the pen- 
tagonal lattice and consider the A^-state clock model on 
it. A brief explanation of the CTMRG method is pre- 
sented. In Sec. Ill we show numerical results on the 
spontaneous magnetization, the internal energy, and the 
specific heat. We summarize the observed phase transi- 
tion. 



II. CLOCK MODELS ON PENTAGONAL 
LATTICE 

We consider the 2D lattice shown in Fig. [1] which is 
a tessellation of regular pentagons. The lattice is in a 
curved plane with a constant negative scalar curvature. 
Therefore, the Hausdorff dimension of the lattice is in- 
finite. For a technical reason in the CTMRG method, 
we have chosen the lattice with the coordination number 
four fl^. Two geodesies drawn by the thick arcs cross 
one another at a site labeled by 9^ . By these two arcs the 
whole lattice is divided into four equivalent parts called 
the quadrants or the corners. 
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FIG. 1: The pentagonal lattice drawn in the Poincare disc. 
The open circles represent the A^-state spin variables . Two 
geodesies drawn by thick arcs divide the system into four 
equivalent quadrants. 



Let us introduce the A''-state clock model on the pen- 
tagonal lattice. On each lattice site there is an A^-state 
spin variable 9^ where i is the site index. The possible 
values of 9.^ are 2tt^/N with ^ = 0, 1, 2, . . . , iV - 1. We 
consider the angle 9^ as the internal degree of freedom. 
Therefore, 9^ has nothing to do with the lattice geome- 
try. If there are only ferromagnetic interactions between 
neighboring spin pairs, the Hamiltonian of the A'^-state 
clock model is written as 



H — — cos I 



(1) 



where J > is the coupling constant. The summation 
runs over all the nearest-neighbor pairs (ij). The case 
A'^ = 2 is nothing but the Ising model with coupling in- 
teraction J and this case has been studied [H, [l^ . The 
case = 4 can be reduced to the Ising model with the 
couphng J/2. We thus chiefly discuss the case A^ = 3, 
which is equivalent to the 3-state Potts model, and the 
cases A^ > 5 in the following. In order to observe the 
phase transition deep inside the system, we impose the 
ferromagnetic boundary conditions so that all the spin 
variables at the system boundary are aligned in the di- 
rection 9 = 0. 

For convenience we represent this clock model as a spe- 
cial case of the interaction-round-a-face (IRF) model on 
the hyperbolic lattice. For instance, let us label the spins 
around a pentagon as shown in Fig. [TJ The IRF weight 
W, which is the local Boltzmann weight corresponding 
to this pentagon, is obtained as 



Wi9, 9,9,9^9,) =Y[exp- 



Jcos (9,-9,+,) 
2fcBT 



, (2) 



where 9q = 9i. Having the IRF weight W thus defined, 
we can express the partition function of the whole system 



{6} 



(3) 



where the product is taken for all the IRF weights in the 
pentagonal lattice. The sum J2{e} taken over all spin 
configurations. 

In order to discuss the phase transition on the hy- 
perbolic lattice, let us consider a system whose size (or 
diameter) L is far larger than the correlation length ^. 
We divide the system into two parts, the boundary area 
(BA) and the deep inside area (DIA). The former, BA, 
is a ring-shaped area, where all the sites in the area are 
within the distance of the order of ^ from the system 
boundary. The latter, DIA, is the rest of the system, 
which we analyze in the following. Because of the hyper- 
bolic geometry, the portion of the BA with respect to the 
whole system is always finite even in the limit L — > oo. 
The situation is similar to that of the Cayley tree [§]. 
Thus the thermodynamic property of the whole system 
is always affected by the boundary condition, especially 
in low temperature l8l|. When ^ is finite, it is possible 
to consider the thermodynamics of the DIA, discarding 
the thermodynamic contribution from the BA, since we 
have assumed L ^ ^ and therefore the size of the DIA 
is sufficiently large. When we collect numerical data of 
the DIA, we always treat sufficiently large systems that 
satisfy L ^ choosing such temperatures for which ^ is 
at most of the order of 1000. We then detect the phase 
transition in the DIA by extrapolation from both low- 
and high-temperature sides. 

We introduce Baxter's corner transfer matrix (CTM) 
C, which represents the Boltzmann weight of a quadrant 
of the system [9]. The partition function Z is then ex- 
pressed as Tr C^, i.e., as the trace of the density matrix 
p = . Applying the concept of the density matrix 
renormalization 2y, [2ll . [2^ , a precise approximation of 
Z can be obtained for large scale systems by way of it- 
erative numerical calculations [l^, [la, ■ These are the 
outline of the CTMRG method, which can be applied to 
statistical models on hyperbolic lattices [3 [3 ■ 

After we obtain the density matrix p for a sufficiently 
large system, we can calculate the expectation values at 
the center of the system, which represent the thermo- 
dynamics deep inside the system. For example, we can 
obtain the spontaneous magnetization 



TWW =Tr [cos(6lJp] /Trp, 



(4) 



where 9^ represents the spin at the center of the system, 
and the internal energy per bond 



-JTr [cos{9^-9'^)p] /Trp, 



(5) 



where 9'^ is the neighboring spin next to 9^. The specific 
heat C*-^' can be obtained by taking the numerical deriva- 
tive of S'^^^ with respect to temperature T. It should be 
noted that J^'^^\ E''^\ and are not thermodynamic 
functions of the whole system but are those of the area 
deep inside the system. 

It has been known that the decay of the density ma- 
trix eigenvalues is very fast for models on the hyperbolic 
lattices [l^, [l^. The clock model under study has the 
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TABLE I: The transition temperatures Tq , the critical ex- 
ponents /3, and positions of the specific heat maximum . 
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1.3659 
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1.3625 


0.5 


0.62948 
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1.3623 


0.5 


0.46295 
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1.3622 
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0.35676 
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1.3622 
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0.28357 
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1.3622 
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0.22997 


13 


1.3622 
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0.13761 
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0.05864 
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1.3622 


0.5 


0.02600 



FIG. 2: Temperature dependence of the spontaneous magne- 
tization M^'^'> for 3 < iV < 30. The open circle denotes the 
discontinuity in Ai^-^\ 



same feature in common. Therefore, it is sufficient to 
keep a very small number of the degree of freedom for 
the block spin variable m in the formalism of CTMRG. 
Typically, we keep m « 2N states. We checked that fur- 
ther increase of m does not improve numerical precision 
in A^(^) and any more, even at the vicinity of the 
phase transition. 



III. NUMERICAL RESULTS 

Throughout this section, we take the coupling constant 
J in Eq. ([T]) to the unit of energy. For all the cases 
N > 2, we observe phase transition, where the transition 
temperatures Tq^' are listed in Tabic HI Note that Tg^' 
converges to Tq°°^ very fast with respect to N. 

Figure [2] shows the spontaneous magnetization Al^^) 
with respect to the rescaled temperature T/T^^\ (Under 
this rescaling, and M^^^ are identical.) If iV = 3, 

the magnetization is discontinuous at Tg . The 3-state 
clock model, which is equivalent to the 3-state Potts 
model, exhibits the first order phase transition if the sys- 
tem is on the pentagonal lattice. This is a kind of mean- 
field behavior, since it is well known that the mean-field 
approximation applied to the 3-state Potts model on 2D 
lattices show the first order phase transition [2^. In the 

vicinity of Tg^' the magnetization A^(^) rapidly con- 
verges to the large N limit M^^\ The inset of Fig. [2] 
displays the low-temperature behavior of M.^'^^ in de- 
tails. Note that in the limit N ^ oo the magnetization 
A^(^' decreases linearly with T at very low tempera- 
tures. Figure [3] shows the square of Ai^^"^ with respect 
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FIG. 3: Square of MS^'^ with respect to (T^^^ - T)/T^ 



to t = (Tg^^^ -T) /Tg^^^ for the cases other than = 3. It 
is obvious that the scaling relation M.'^^'^ oc is satisfied 
with the exponent /? = i. 

Figure [4] shows the internal energy f There is a fi- 
nite jump in E^'^^ at Tg'^^ where the latent heat per bond 
C — — fi^"* is 0.078. Analogously to the magneti- 
zation M.'^^\ the f is linear in T at low-temperature 
region in the limit TV — > oo. 

Figure [5] shows the rescaled specific heat C^^^/C^mlx, 



where C^max is the specific heat at T^^^ , with respect to 
the rescaled temperature T/Tq^\ Evidently, a disconti- 
nuity in the specific heat is observed for the cases N = 2 
and TV > 4. Thus, the second order phase transition 
has the mean-field nature. There is no indication of the 
BKT transition that is observed for clock models on flat 
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FIG. 4: The absolute value of the internal energy l^'^']. The 
open circles denote the jump in the case N = 3. 



0.7 
0.6 
0.5 
0.4 
'o.3 
0.2 
0.1 
OJD, 



PN=13 



Q'N=2Q 

'orooo 0.005 0.010 0.015 0.020 0.025 o.oso 

FIG. 6: The Schottky peak position Tg'^' versus 



1.0 

0.9 

0.8 

0.7 

g gO.6 
U 

^ 0.5 
U 0.4 
0.3 
0.2 



1.5 

, 1.0 



0.5 



0.0 



.1 1 1 1 1 1 1 1 1 1 1 1 


C(C): 







0.0 



_30 13 



0.5 



1.0 
T 



-120/ 10 / 



ill,'//.'' / 

Q J. Li-it: Cii ' 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 1.1 



T/T. 



(N) 



FIG. 5: The rescaled specific heat C'^^'/C' 

max versus the 

rescaled temperature T/T^^K The inset shows a typical ex- 
ample for the case A*' = 10 without rescaling. 



2D lattices [2||. 

When N is larger than 5, we observe the Schottky 
type peak in the specific heat. Figure [6] shows the N 
dependence of the Schottky peak position Tg^-*. As it 

is shown, Tg^^ is proportional to 1/N'^. This is qualita- 
tively in accordance with the energy scale of local excita- 
tion 2{2'K / NY J from the completely ordered state. It is 
thus concluded that the Schottky peak disappears in the 
limit iV oo and that the specific heat of the classical 
XY model on the pentagonal lattice remains finite even 
at T = 0. 



IV. CONCLUSIONS 



We have studied the A^-state clock models on the pen- 
tagonal lattice, which is a typical example of the hyper- 
bolic lattices. The phase transition deep inside the sys- 
tem is observed by use of the CTMRG method. From the 
critical exponent /3 = ^ for the spontaneous magnetiza- 
tion and the jump in the specific heat, we conclude that 
the phase transition for = 2 and > 4 is mean-field 
like, provided that the ferromagnetic boundary condi- 
tions arc imposed. The Hausdorff dimension, which is 
infinite for the hyperbolic lattices, is essential in the ob- 
served critical behavior. We conjecture that the phase 
transition deep inside the system is also present for sys- 
tems with free boundary conditions. 

In the case when A^ = 3, where the system is equivalent 
to the 3-state Potts model, we observed the first-order 
phase transition. Since the q-state Potts model tends 
to exhibit the first-order transition for larger q [2^, it 
is expected that the transition of g > 3 Potts models 
on the pentagonal lattice is of the first order. We have 
partially confirmed the behavior for several values of q 
and we conjecture that the transition is of the first order 
on any kind of hyperbolic lattices when (7 > 3. 

We observed stable ferromagnetic states below Tjj^'' 
even in the continuous limit A" — > 00. This fact does 
not contradict to the Mermin- Wagner theorem [l| since 
the pentagonal lattice is not on the flat 2D plane. The 
vortex energy on hyperbolic lattices might be larger than 
that on the flat lattice. The difference may elucidate the 
absence of the BKT phase transition on the pentagonal 
lattice. 
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